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Abstract. The ordered alloy La2Lio.5Cuo.5O4 is found to be a band insulator in local 
density approximation (LDA) calculations with the unoccupied conduction band having 
predominantly d x i _ y 2 -symmetry and substantial weight in O 2p-orbitals. This is equiva- 
lent to a predominant local singlet configuration d 9 L or a low spin Cu 3+ -ion with both 
holes in orbits having d x 2_ y 2 -symmetry, i.e. Zhang-Rice singlets. A fairly modest reduc- 
tion of the apical Cu-0 bondlength is sufficient to stabilize a high spin triplet Cu 3+ -ionic 
configuration with holes in both d x 2_ y 2 and d 3z 2_ r 2 orbits in LDA + U calculations. This 
leads us to identify the low energy triplet excitation found in NQR studies by Yoshinari et 
al. as a local high spin Cu 3+ -ionic configuration accompanied by a substantial reduction 
of the apical Cu-0 separation, i.e. a anti-Jahn- Teller triplet polaron. 
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The substitution of lithium for copper in La2Cu04 is formally equivalent to Sr doping 
on the La site as each donates one hole per dopant. In lightly doped samples Li and Sr 
doping have essentially indistinguishable effects on the magnetic properties and lattice 
structure |],|] . However in-plane substitutions for Cu and out-of-plane substitutions 
for La were found to be very different in the conductivity. Li +1 has ionic radius essen- 
tially the same as that of Cu +2 , and brings a hole with it into the plane. But the alloys 
La2Cui_ a; Lij 3 ;04 are never metallic nor superconducting so that this hole must be localized 
in contrast to the mobile holes introduced by out-of-plane substitutions. The solid solu- 
bility of Li is such that complete filling of the copper band ( at 50% Li) can be achieved. 
At this composition the Li and Cu ions form an ordered superlattice || in which all 
Cu ions are surrounded by four in-plane Li +1 ions (Is 2 , closed shell electronic configura- 
tion), leading to isolated Cu0 4 clusters. This compound was found to be a diamagnetic 
insulator. 

The Nuclear Quadrupole Resonance (NQR) study of La2Cuo.5Lio.5O4 [|J reveals a 
magnetic excitation of the doped-hole state with an energy of ~ 130 meV. This energy is 
much smaller than the estimation of the singlet-triplet splitting of the local Zhang-Rice 
singlet m which is of the order of few eV. There are also indications that this magnetic 
excitation is coupled to the charge fluctuation or lattice distortion around the Cu site 0. 

La2Cuo.5Lio.5O4 is not the only compound with formally trivalent copper which is a 
diamagnetic insulator. Another example is NaCu02 @. Spectroscopy measurements and 
Configuration Interaction (CI) calculations by Mizokawa et al. |7j led them to conclude 
that the ground state of NaCu02 is dominated by d 9 L (L: ligand hole) configurations and 
not by d 8 as for a simple Cu 3+ ion. The ligand holes L are not metallic because strong 
p — d hybridization between d 8 and d 9 L configurations leads to a split-off rf 9 L-like discrete 
state above the oxygen continuum. 

Calculation of the electronic structure of NaCu02 in Local (Spin) Density Approxima- 
tion (L(S)DA) || showed that this compound can be described as a conventional insulator 
with a band gap arising from simple covalent effects {p — d hybridization). The unoccu- 
pied band has more oxygen character than copper thus confirming the conclusion from CI 
calculations that "doped" holes are situated predominantly on oxygen orbitals. The key 
structural element of the crystal structure of NaCu0 2 is the same as in layered cuprates 
(for example La2Cu04): CUO4 plaquettes with the Cu atoms in the center of a square 
of oxygen atoms. However in layered compounds the squares are corner-sharing and the 
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Cu-O-Cu bond angles are nearly 180° resulting in very broad pda-band, while in NaCu02 
the squares are edge-sharing and those angles are close to 90° and the corresponding 
bandwidth is much smaller. 

In La2Cuo.5Lio.5O4 the CUO4 plaquettes are separated by Li ions, so that there is no 
Cu-O-Cu bonds at all and one can expect an even narrower pda-b&nd than in NaCuCV 
We have performed LDA calculation for La2Cuo.5Lio.5O4 using the linearized muffin-tin 
orbital method in an atomic-spheres approximation (LMTO-ASA) ||. The results (Fig.l) 
show that, indeed, the ground state is a nonmagnetic insulator with a sizable gap value 
(« I eV), compared to the 0.3 eV value in NaCu0 2 0. The unoccupied band is again 
rather narrow and with a symmetry of Cu centered d x 2_ y 2 orbitals. It contains 40% Cu3d 
orbitals and 60% 02p states, indicating the strongly covalent nature of the singlet ground 
state. The top of the valence band is predominantly oxygen in origin, however there is 
significant admixture of d^ z 2_ r 2 orbitals of copper. 

Our results shows that in order to reproduce singlet ground state of the doped hole in 
La2Cuo.5Lio.5O4 there is no need to take into account Coulomb interaction corrections to 
the one-electron LDA. However the magnetic excited state cannot be explored in a calcu- 
lation scheme which does not include Coulomb interactions inside the <i-shell of copper. 
For example the undoped cuprates experimentally are antiferromagnetic insulators while 
LDA gives a paramagnetic metallic ground state. The reason for this discrepancy is that 
while in LSDA the splitting between majority and minority spin states is driven only by 
the exchange interaction with a value of the Stoner parameter of ~ 1 eV, the real driving 
force for the antiferromagnetic insulator solution must be the much larger direct Coulomb 
interaction parameter U ~ 8e V. This contradiction was resolved in the so-called LDA+U 
method where orbital-spin polarization caused by the Coulomb interaction is directly 



taken into account 10,11 



The main idea of the LDA + U method is that LDA gives a good approximation for 
the average Coulomb energy of d-d interactions E av as a function of the total number of 
(^-electrons N = J2 ma n ma where n ma is the occupancy of a particular d mtr -orbital: 

E av = \U N(N - 1) — I J N(N - 2). (1) 

But LDA does not properly describe the full Coulomb and exchange interactions between 



(^-electrons in the same d-shell. So Anisimov et al. |10 , |TT| suggested to subtract E av from 



the LDA total energy functional and to add orbital- and spin-dependent contributions to 
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obtain the exact (within a mean-field approximation) formula: 

E = E LDA — E av + | U, 



mm' ^mcr^m' ' —a 

m,m' ,cr 



+ 2 ^ ^ (Umm' J mm'^^ma^m' a ■ (2) 
m^m' ,m' ,cr 

Taking the derivative w.r.t. n ma gives the orbital-dependent one-electron potential: 

V m a{r) = V LDA {r) + J](l7 mm < - U eff )n m/ _ a 

m! 

+ ^ {Umm 1 ~ Jmm 1 ~~ U e ff)n ma (3) 
m'j^m 

+ Ueffiz ~ n "ia) — \J ■ 

with U e ff = U — \ J ■ The Coulomb and exchange matrices U mm i and J mm i are expressed 
through the integrals over products of three spherical harmonics and screened Coulomb 
and exchange parameters U and J |fL0| . 

A nontrivial problem is what value of the screened Coulomb interaction U to use. For 
insulators, such as late transition-metal oxides a good approximation is to calculate static 
screening of the d-d intrashell Coulomb interaction in a supercell LDA calculation |T2 . 



The question is what symmetry should the lowest energy excited states have? The 
band gap separates states which are both mainly oxygen but due to the hybridization with 
Cu3d orbitals the states have symmetry of x 2 — y 2 and 3z 2 — r 2 for the unoccupied and 
occupied bands respectively. The x 2 — y 2 band is higher in energy due to the Jahn- Teller- 
distorted CuOe octahedra in La2Cuo.5Lio.5O4 crystal structure: the length of the Cu-0 
bond in the ab-plane is 1.8 A but the distance to the apical oxygen is larger - 2.4 A B. It 
follows that the lowest energy excitation will be from the (formally) d x 2_ y 2^d x 2_ y 2\ i singlet 
configuration to the d x 2_ y 2^d^ z 2_ r 2^ triplet configuration. 

We have performed LDA+U method calculations for both singlet and triplet config- 
urations. We find that a starting triplet configuration is not a stable solution but self- 
consistently converges to the singlet solution. (The singlet solution of LDA+U method is 
practically the same as for pure LDA, because in the absence of the orbital-spin polariza- 
tion the LDA+U correction to LDA is irrelevant.) However, if the copper- apical oxygen 
bond length is contracted by 16% (0.38A), then a stable triplet solution appears and be- 
comes the ground state. For 16% contraction the total energy of the magnetic solution is 
still 0.06 eV higher than the nonmagnetic one, however already for 17% contraction the 
triplet energy is 0.1 eV lower. 
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It is instructive to follow the changes, with copper-apical oxygen bond contraction, 
of the two bands below and above Fermi energy (Fig. 2). For the undistorted structure 
(Fig. 2a) those bands (the lower one of d 3z 2_ r 2 symmetry and the higher one of d x 2_ y i) are 
well separated from each other and the admixture of the c^.^-orbitals to the valence 
band is relatively small. With distortion the hybridization of the Cu3c?3 2 2_ r 2-orbitals with 
the apical oxygen 2p 2 orbitals becomes stronger and the antibonding band goes higher in 
energy and the admixture of the dj, z i_ r i -orbitals in this band becomes stronger (Fig.2b). 
In order for the magnetic state to be stable, the splitting between spin-up and spin- 
down bands must become large enough that the top of the spin-up d x 2_ y 2 band does not 
overlap with the bottom of spin-down d? jZ 2_ r 2 band (in other words magnetic energy must 
overcome kinetic energy). The Fig.2c shows that for the critical value of the distortion 
(16%) this condition is nearly satisfied. 

Similar effects were found in a supercell LDA+U calculation for a doped hole in 



La2- x Sr x Cu04 ||13|| . In that work two solutions were found for a hole introduced in the 
CuC>2 plane: one solution had the symmetry x 2 — y 2 with the hole spin antiparallel to 
the d-hole of Cu atom, while the other one had 3z 2 — r 2 symmetry and spin parallel to 
the Cu spin. The latter solution was present only for a contracted apical Cu-0 distance 
(0.26Acontraction). The total energy calculation in a full-potential scheme including the 
lattice relaxation showed that the total energy minima for these two solutions were very 
close in energy with the triplet state only 54 meV higher than the ground state singlet. 

As the hopping between CuOg octahedra in La2Cuo.5Lio.5O4 is smaller than in 
La2Cu04, one would expect more localized states in the former and, hence, a larger 
separation between the singlet ground state and excited triplet state. The value found in 
NQR measurements j|] 130 meV looks quite reasonable from this point of view. 

We would like to emphasize that while the excited magnetic state is a triplet with 
S=l, the actual magnetic moment residing in Cu3a? orbitals found in LDA+U calculation 
is very small: 0.8 That is due to the fact that the bands (or Wannier orbitals) of 
x 2 — y 2 and 3z 2 — r 2 symmetry have only 40% of Cu3d orbitals contribution, being mainly 
oxygen in origin as implied by the predominance of d 9 L many electron configurations in 
model CI calculations 0. 

In conclusion, we find that LDA correctly gives a band insulator ground state for 
La2Lio.5Cuo.5O4 with the unoccupied conduction band having predominantly d x 2_ y 2- 
symmetry and substantial weight in O 2p-orbitals - a state equivalent to a formal va- 



5 



lence Cu 3+ in a low spin configuration or a Zhang- Rice singlet. A reduction of the Cu-0 
apical distance stabilizes a high spin configuration and this leads us to identify the 130 
meV triplet excitation observed in NQR with such a local configuration i.e. an anti-Jahn- 
Teller polaron. Two aspects of the experiments remain to be clarified - the origin of 
the low energy nuclear spin relaxation process which dominates at lower temperature at 
T < 170K and secondly the absence of a significant activated contribution to the uniform 
susceptibility from the 130 meV triplet excitations - at least below room temperature. 

The work was partly supported by Russian Basic Research Foundation (RFBR grant 
96-02-16167). One of us (VIA) wishes to thank the 'Zentrum fiir Theoretische Studien' 
at the Institute for Theoretical Physics, ETH-Ziirich for hospitality. 
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FIGURES 

FIG. 1. The total (a) and partial (02p (b) and Cu 3d (c)) densities of states (DOS) for 
La2Lio.5Cuo.5O4 calculated in the undistorted crystal structure. The bands formed by different 
d-orbitals are marked by arrows. 

FIG. 2. The partial Cu 3d densities of states (DOS) for La2Lio.5Cuo.5O4. Upper panel (a): 
nonmagnetic solution in the undistorted crystal structure. Middle panel (b): nonmagnetic solu- 
tion with 16% contraction of the copper-apical oxygen bond length. Lower panel (c): magnetic 
solution with 16% contraction of the copper-apical oxygen bond length. 
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